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ABSTRACT 

We perform a time-resolved spectral analysis of GRB 130606B within the framework of a fast-cooling 
synchrotron radiation model with magnetic field strength in the emission region decaying with time, 
as proposed by Uhm & Zhang. The data from all time intervals can be successfully fit by the model. 

The same data can be equally well fit by the empirical Band function with typical parameter values. 

Our results, which involve only minimal physical assumptions, offer one natural solution to the origin 
of the observed GRB spectra and imply that, at least some, if not all, Band-like GRB spectra with 
typical Band parameter values can indeed be explained by synchrotron radiation. 


1. INTRODUCTION 

Despite the multifarious observations over 40 years 
since their discovery, the prompt emission phase of 
Gamma-Ray Bursts (GRBs) is still not fully understood 
in several fundamental aspects (e.g., |Zhang| |201l] for 
a review): (1) What is the composition of the ejecta 
(baryonic fireball or Poynting flux dominated jet)? (2) 
What are the energy dissipation and particle accelera¬ 
tion mechanisms (from shocks or via magnetic recon¬ 
nections)? (3) What is the radiation mechanism (syn¬ 
chrotron, synchrotron self-Compton, or Comptonization 
of quasi-thermal emission from the photosphere)? Inter¬ 
pretation of the observed GRB energy spectra plays a 
key role to answer these questions. A direct identifica¬ 
tion of the GRB prompt emission mechanism would lend 
further constraints on the energy dissipation mechanism, 
and eventually on the composition of the GRB jets. 

GRB spectra have traditionally been described by a 
joint smoothly broken powe r law function - the “Band” 


function (Band et al. 1993). This function successfully 


Briggs et al. 

1999; Abdo et al. 2009 Zhang et al. 2011) 

but it has been difficult to associate i 
the physical processes responsible foi 
Two leading candidate radiation m< 
prompt emission are synchrotron i 

the parameters with 
: the GRB emission, 
echanisms for GRB 

'adiation (Meszaros 

et al. 

1994 

Daigne & Mochkovitch 

1998| Lloyd & Pet- 

rosian 

2000 

Zhang & Yan 2011; I 

Jaigne et al. 2011) 

and Comptc 

photosphere 

inization of quasi-thermal photons from the 

of the fireball (ffhomi 

pson 1994 Rees & 

Meszaros 20051 Beloborodov 2010; I 

jazzati & Pegelman 


values of E p (peak energy in the vF v spectrum) and 
P (the high-energy power-law spectral index). The 
straightforward predictions of the low-energy power-law 
index a of both models, however, deviate from the typical 
value a ~ — 1 revealed by the observations (e.g., |Preece 


et al. 2000 Zhang et al. 2011 


Nava et al. 


2011). The 


traditional fas t-cooling synchro tron model predicts an a 
value of —1.5 ( Sari et ah] [l998). This has been regarded 


as a strong drawback of the traditional fast-coolin g syn¬ 
chrotron model t o interpret GRB prompt emission (Ghis- 
ellini et al.|2000). The photosphere models, on the other 
hand7predict too hard on spectrum. A pure blackbody 
emission in Rayleigh-Jeans regime predicts a ~ +1. Rel¬ 
ativistic effects (equal-arriv al-time surface effect) would 


soften it to a 

|2 

E 


+0.4 (e.g., Beloborodov 2010 Deng & 


Zhang 12014), which is still too hard compared to the 


data. Effects that modify the low-energy photon index 
closer to -1 include adding synchrotron self-absorption 
or a smoot h cutoff in the Jow-energy electron distribu¬ 
tion (Lloyd & Petrosian 2000), addin g a ph otosphere 
comp onent in the synchrotron model (Meszaros & Rees 


sphere model ( 

Toma et al. 

2011; Vurm et al. 

20n|), or 

introducing a certain type ol structured jet 

Lundman 


et al. 2013). However, some of these modifications ei¬ 
ther invoke simplified assumptions or introduce spectral 
features that are inconsistent with the observations. 

Observations of GRBs with the Fermi Gamma-Ray 
Observatory have revealed additional features beyond 
the B and function, including a quasi-therm a l compo- 
nent (|Ryde et al. 2010 Guiriec et al. 2011 |Axelsson 


et al. |2012[ Guiriec et al. 2015[ ) and a power-law com- 


Zhang et al. 

2011j. Through a detailed time-resolved 

spectral analysis, Zhang et al. 

(2011) suggested that the 


Band component, the quasi-thermal component, and the 
high-energy component are three elemental spectral com¬ 
ponents in GRB spectra. A recent d etailed stud y of th e 
time-resolved spectra of three GRBs (Guir iec et al.|2 015) 
indeed suggests that the interplay among three elemen- 
tal components can explain a variety of observed GRB 
spectra. These observations have spurred efforts to ap¬ 
ply physical modeling to the GRB data and uncover the 
physical meaning of various spectral components. The 
quasi-thermal component is often interpreted as em issio n 
from the photosphere of the outflow (Meszaros & Rees 
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2000). Since this component is predicted to be bright in 
the standard fireball models, non-detection or a detec¬ 
tion of a weak photosphere co mponent point toward s a 
Poynting-flux-dominated flow (Zhang & Pe’er 2009) or 













































































































































Table 1. Sample and Fitting Results 
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a hybrid outflow with significant Poynting flux in the 
jet (Gao & Zhang 2015). With the quasi-thermal com¬ 
ponentIdentified as the* photosphere emission, the main 
Band component is more l ikely of a synchrotron origi: 


Burgess et al. I (2011 2014) found that the addition ofa 


photospheric black body to a synchrotron radiation com¬ 
ponent reproduces the observed spectrum of eight bright 
GRBs detected by Fermi , provided that the electrons are 
in the slow cooling regime. This is however contrary to 
expected, since in the GRB emission region, electrons 
likely cool rapidly due to a strong local magnetic field. 
The simplest fast-cooling spectra, even with the super¬ 
position of the photosphere component, were found in¬ 
compatible with the observations. 

In these and previous studies, the magnetic field in the 
GRB emission region is assumed to be constant. How¬ 
ever, in a rapidly expanding jet source such as a GRB, 
flux conservation indicates that its magnetic field can¬ 
not be preserved as a constant. This important physi 


cal co nsi derat ion has been taken into account by Uhm & 
Zhang p014). By introducing a decrease in the magnetic 


held witlTincreasing radius as a (conical) jet expands into 
space and by assuming that the GRB emission region is 
far from the central engine (typically 10 15 cm and be¬ 


yond), they create the following picture (Uhm & Zhang 
2014): the GRB emission is assumed to be generated 


from fast-cooling synchrotron radiation of electrons ac¬ 
celerated in the emission region (either internal shocks or 
internal magnetic reconnection sites), which is relatively 
far away from the central engine. Electrons are continu¬ 
ously injected as a power law distribution (N(^y) oc y _p ) 
with a minimum Lorentz factor 7 i n j, and with a certain 
injection rate as a function of time Ri n j = R^ n -(t/to) q . 
Since the emission region is streaming outwards in an 
expanding jet, the co-moving magnetic held strength de¬ 
cays with radius as B'(r) = B' 0 (r/ro)~ b . The accelerated 
electrons undergo synchrotron cooling in the co-moving 
frame. Since the cooling time scale is shorter than the 
dynamical time scale, the electrons are in the fast cool¬ 
ing regime, so that an electron population below 7 i n j is 
developed as a function of time. In the standard fast 
cooling model where B' is assumed a constant, the elec¬ 
trons below 7 i n j has a spectrum of N( 7 ) oc 7 “^ with 
p = 2, which gives a photon index a = —1.5, inconsis¬ 
tent with the data. However, in a decaying B' held, the 
degree of cooling varies as a function of time, so that the 
electron spectrum below 7 i n j is curved with an effective 
p < 2. The photon spectrum therefore hardens. Varying 
the parameters, the predicted effective a can be as hard 
as - 0 . 8 , with the typical value a ~ — 1 reproduced for a 
range of parameters. This concept for the hrst time in¬ 
troduces a physical fast-cooling synchrotron model that 
accounts for the typical observed a value^] 

The investigation of this concept has thus far been lim- 


2 Another support to this view is from GRB 110721A, which 
showed an extremely high E p 
the temperature of the photosphere 1 
2012 ). 


view 10 iiu 1 iiir 

(|Axelsson et al. 2012 | far exceeding. 
phere (Zhang et al. 2012 Veres et al. 

’ Zhao et al.] ([2014) proposed a different scenario by invoking 
decay of magnetic held locally behind the shock front within the 
framework of the internal shock model. They showed that by con¬ 
sidering both synchrotron and inverse Compton cooling of elec¬ 
trons, the typical a ~ — 1 can also be reproduced within that 
model. With a small emission radius, this model may struggle to 
provide a natural explanation for the seconds-long pulses observed 
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ited to numerical calculations that have produced spec¬ 
tral parameters and temporal evo lution properties sim - 
ilar to those observed in GRBs (Uhm & Zhang 2014), 
but the model has not been applied to real data, in this 
paper, we aim to test whether such a fast-cooling syn¬ 
chrotron model can really account for the observed GRB 
spectra. We describe our data reduction method in §2 
and fitting tools in §3. Our fitting results are presented 
in §4, followed by a summary and discussion in §5. 

2. DATA REDUCTION 


As will be discussed in £4.2, it is computationally ex¬ 


pensive to fit the observed spectra with the numerically 
calculated synchrotron model. We select one typical 
GRB detected by Fermi/ GBM as a test, whose time- 
resolved spectra can be well fit by the Band function. 
Because most GRBs have strong s pectral evolution dur¬ 
ing the prompt emission phase (e.g. Zhang et al.|2011 Lu 


p nase y 

et al.|2012 Gniriec et al.||2015| ), this work focuses on the 
time-resolved data. We select the second brightestnGRB 
in the second GBM catalog, i.e., GRB 130606B (GBM 
trigger number 392212536/130606497). The brightness 
of the GRB 130606B enables us to perform a detailed 
time-resolved spectral analysis with fine time-resolution, 
owing to the large number of y-ray photons. This GRB 
was about 91 deg from the LAT boresight at the time of 
the trigger and remained well outside the nominal field 
of view for the standard data analysis during the whole 
prompt emission 100 s). A LAT detection during the 


(LLE) data ( 

Vianello et al. 

2013 

), but the energy re- 

sponse of the LLE data at tl 

ris large offset angle to the 


boresight is not reliable. We thus omit LAT data from 
our analysis. 

The Fermi/ GBM data are pr ocessed usi ng the stan¬ 
dard procedure described in Zhang et al.| (2011). We 
selected the data from the two brightest Nal detectors 
and the brightest BGO detector. The light curve of the 
brightest Nal detector (Nal 8) is shown in the bottom 
panel of Figures [l] and [2] Our analysis only focuses on the 
main prompt emission episode (from 0 s to 59 s; roughly 
the episode covered by Tgo) of this burst. Based on the 
signal-to-noise ratio, we divide this time interval into 24 
slices (Table 1 and Figures |T] and [2]). For each slice we 
extract the source spectrum, background spectrum and 
the corresponding response files, which are needed to per¬ 
form the spectral fitting. 


3. A GENERAL MONTE-CARLO FITTING APPROACH 

Generally speaking, performing a fit does not only 
mean finding the “best-fit” parameters, that bring a par¬ 
ticular model into optimum consistency with the data, 
but also means, more importantly, to find the true like¬ 
lihood map of the entire parameter space, with which 
one can know how the model can be really constrained 
in a global confidence manner. The latter part is more 
challenging since theoretical models always involve multi¬ 
ple model parameters that may be subject to correlation 
and multi-modality of the likelihood in their parameter 
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Figure 1 . Bottom: Light curve of GRB 130606B in 15-350 keV of 
Fermi GBM detector Nal 8. Red numbers mark the slice intervals 
listed in Table 1. Top 1-3 panels : Evolution of the Band function 
parameters via spectral fits of each time interval. 


space. Traditional fitting algorithms (e.g., Levenberg- 
Marquardt) typically fail to map the multi-dimensional 
likelihood. We have developed a fitting engine ( McFit ) 
that employs a Bayesian Monte-Carlo (MC) fitting tech¬ 
nique to reliably fit parameters that are constrained by 
the data even when other parameters are unconstrained. 
Such a technique has been successfully applied t o Swift 
afterglow modeling of GRBs with 7 parameters (Zhang 
et ah| 2015). Using this technique, the best-fit parameters 


and their uncertainties can be realistically determined by 
the converged MC chains. 

To compare a theoretical spectral model to data, we 
need to convolve the model with instrumental response, 
the Detector Response Matrix (DRM). We have written 
a general-purpose forward-folding code that can handle 
the instrument response R(J, E ) from any mission (e.g. 
Fermi/ GBM), read in any model spectrum, F(E,P), re¬ 
gardless of its form (e.g. a mathematical function or 
a tabulated array), and calculate the model-predicted 
count spectrum Cm(G P) = F(E, P)R(I, E)dE, 
which can be directly compared with the observed count 
spectrum C(I). We then calculate the likelihoodj for 
those Cm(I,P) and C(I) pairs. Such a likelihood func¬ 
tion is then input into our fitting engine McFit so that 
the parameter space can be explored. A final software 
package, McSpecFit , combines this Bayesian MC engine, 
the general forward-folding codes and the likelihood cal¬ 
culations, and can fit any models to the observed spectra. 


4. RESULTS 

4.1. Band Function Fits 


in GRBs as well as the hard-to-soft and tracking behaviors of E p 
within these pulses. 

4 The brightest one is 130427A, which is subject to saturation 
in the GBM TTE data, and hence not selected. 


5 In this paper, we use the maximum likelihood-based statistics, 
the socalled PGSTAT, for Poisson data, given by Cash (1979), 
but assuming the background to be Gaussian; see also https:// 
heasarc.gsfc.nasa.gov/xanadu/xspec/ 
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We first fit the observed time -resolve d spectra with the 
empirical Band function (Band Tet al.|[l9 93): 


F(E) = A 


where 


( lOO E kev )°^ E<(oc-P)E 0 , 

{toAsvY ( ^Oollv 0 )° P ’ 

( 1 ) 

= (2 + ol)Eq, ( 2 ) 


A is the normalization factor at 100 keV in units of pho¬ 
tons s -1 cm -2 keV -1 , a is the low-energy photon spec¬ 
tral index, /3 is the high-energy photon spectral index, 
and E p is the peak energy in units of keV in the observed 
vF v spectrum. 

The results of our fits for each slice are shown in Table 
1. The measurement of the goodness-of-fit in each slice 
is plausible, with PGSTAT/dof ^ 1, as shown in the 
5th column of Table 1. As an example, we show the 
parameter constraints resulting from our McSpecFit in 
the inner plot of Figure [3j We also plot the evolution of 
the three characteristic parameters of the Band function 
model in Figure 1. We notice that except for 9 slices 
(intervals 0, 2, 3, 4, 5, 17 and 21) the values of /3 are 
constrained to have fairly large values (2.5-8) and large 
uncertainties, which indicates that a cutoff power law 
model would also give equally good fits. We therefore 
mainly focus on the properties of a and E p . 

A strong spectral evolution is shown in Figure 1. E p 
and the peak count rate (or equivalently, peak flux) in 
each slice correlate with each other and show a clear 
tracking behavior (see also Ford et al 2015, Lu et al 2012). 


4.2. Synchroton Model Fits 
Next, we apply the fast-cooling synchrotron ^emission 


model with a decaying magnetic field (Uhm & Zhang 


2014) to fit the time-re solved spectra. 


Mathematically, the Uhm & Zhang (2014) model is 
defined by the following nine parameters: 

• the bulk Lorentz factor Y of the emitting region; 

• the magnetic field strength B 0 (in units of Gauss ) 
at 10 15 cm; 

• the magnetic field decay index 5; 

• the power-law index p of the injected electron spec¬ 
trum; 

• the minimum injection Lorentz factor 7 i n j of elec¬ 
trons; 


• the power-law index q of the injection rate as a 
function of time; 

• the normalization of the electron injection rate i2? n j 
in units of s -1 ; 


• the observer-frame time since the first electron in¬ 
jection t in units of s of a given emission episode; 

• the redshift z. 
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Figure 2. Bottom: Same as the bottom panel of Figure 1. Top 
1-8 panels : Parameter evolution of the synchrotron model fitting 
of each slice. 


The theoretical flux density (in unit of mJy) can be 
calculated numerically in the form 

F V {E) = M(E-T,p,y m j,R° ni ,q,B 0 ,b,i,z). (3) 


We im plement the numerical code of |Uhm fc~Z hang 
(2014) into the data fitting code, having it serve as a 
black box to generate numerical model spectra given any 
set of the input parameters. 

The direct output of the modeling is F U (E). It is 
converted to photon flux and then read into McSpecFit. 
Spectral fitting is then performed following the Bayesian 
Monte Carlo fitting method within the allowed parame¬ 
ter ranges (Table 2) as introduced in Section 3, until the 
best-fit model parameters are obtained. 

Two caveats are worth mentioning. First, GRB 
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Figure 3. Parameter constraints of the synchrotron model ( up-right ) and Band function (bottom left) fitting for Interval 7. Histograms 
and contours show the likelihood map of the parameter-constraint outputs by our McSpecFit package. Red crosses mark the best-fit values 
and 1-sigma error bars. All parameters are constrained in reasonable ranges (see Table 1). 


130606B does not have a measured redshift. We assume 
z = 1 in the modeling. The shape of spectrum does not 
depend on the redshift. For a different redshift, the same 
observed spectrum can be reproduced by essentially ad¬ 
justing the injection rate parameter Rf n y Second, in the 

fitting, t is taken as a free parameter. Strictly speaking, 
it should be defined as the observer-frame time since the 
injection of the first group of electrons in a particular 
emission episode. However, since GRB light curves in¬ 
clude overlapped pulses, it is hard to track the number of 
emission episodes. Furthermore, the predicted flux is not 
very sensitive to t if the duration is in the time scale of 


seconds ([uE m & Zhang 2014). This can be readily seen 
in Figj3j which shows that the contour of t is relatively 
broad. We therefore make no effort to guess the t value 
for each time bin, and leave it as a free parameter. 

The results of fitting the data to the synchrotron model 


are presented in Table 1 and Figure [2] Similar to the 
Band function fitting, we show an example of the pa¬ 
rameter constraints for the same time interval 7 in the 
outer plot of Figure [3J We show that the synchrotron 
model provides equally good fits to the data as the em¬ 
pirical Band function model (Figure |4|. 

The evolution of various model parameters is shown 
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Figure 4. Comparison between the PGSTAT of the Band func¬ 
tion fitting and the synchrotron model fitting for the 25 time slices 
in Figure 1. The dashes mark the line of equality. 


Parameter 

Range 

T 

[10 2 , 10 3 ] 

P 

[2,6] 

7inj 

[10 3 , 10 6 ] 

^ nj (a* 1 ) 

[10 44 , 10 50 ] 

q 

[0.01, 5] 

Bo (G) 

[5, 1000] 

b 

[0.9, 2.0] 

t (s) 

[0.1, 10.0] 


Table 2 

Allowed Range of The Parameters of the Synchrotron Model in 
the Monte Carlo fit 

in Figure [2| Similar to the Band function parameters, 
these model parameters all undergo variations through¬ 
out the burst. As expected, the magnetic field strength in 
the emission region is found to decay with radius, with b 
varying in the range between 1 and 2. However, the large 
error of b makes it possible that b does not evolve signifi¬ 
cantly throughout the burst. This is consistent with the 
expectation: whereas the flux and E p vary significantly 
throughout the burst due to the change of electron in¬ 
jection rate and typical Lorentz factors, the geometric 
configuration and the degree of magnetic dissipation of 
the jet (which defines b) do not evolve significantly. 


4.3. Comparison between the Two Model Fits 

In Figure [4] we compare the goodness of fits, PGSTAT, 
of the Band model and the synchrotron model. One 
can see that the two models give essentially identical 
goodness -of-fit. This shows tha t the physical model pro¬ 
posed by |Uhm & Zhang| (|2014|) successfully fits the typi¬ 
cal GRB spectra which are otherwise fitted by the empir¬ 
ical Band function (virtually a mathematical shape). In 
Fig 5j we show a comparison between the photon spectra 
of the two model fits as well as the residuals of the fits for 
the time interval 7 as an example. Again it clearly shows 
that both models are adequate to represent the data. In 


Fig(6j we show the model curves of the Band function and 
synchrotron model. One can see that the two curves are 
esse ntially identical with in the GBM energy band (see 
also Uhm fc Zhang||2Q14 ), and only show noticeable dif- 
ference in a much wider energy band. This explains the 
capability of the synchrotron model to mimic the Band 
function within the GBM band. 

The E p tracking flux feature (Figure [l]) shown in this 


Ford et al. 

1995 

Lu et al. 

2012) 


lution pattern is predicted (Uhm & Zhang 2014|). Our 
analysis suggests that by allowing variation of the pa¬ 
rameters especially 7 i n j, the synchrotron model can also 
reproduce ^-tracking patterns in GRBs. 

The hardest a in our Band function fit is — 0.521^24 
(interval 0). The best-fit value -0.52 already exceeds the 
so-called synchrotron “line-of-death” (Preece et al. 1998). 
However, due the large error of <a, the data are also con- 
sistent with the synchrotron model. This suggests that 
“line-of-death” is not a hard limit for the synchrotron 
model. 


4.4. What makes the Band Spectrum? 

A relevant question is what physical processes shape 
a synchrotron spectrum to be “Band-like”. In order to 
address this point, we systematically checked the pair 
correlations between the fitted Band parameters and the 
synchrotron model parameters. Among those correla¬ 
tions, we find the following two pairs (Figure 7) showing 
“strong” (Dancey & Reidy 2004) correlations with Pear¬ 
son correlation coefficients, \r\ >0.7: <a-7i n j and <a-g; and 
the following four pairs showing “moderate” correlations 
with 0.4 < |r| < 0.6: <a-T, <a-H 0 , E p - 7mj, E p -R® nj . This 
indicates that the observed spectral properties can be 
mostly accounted for by the variations of the jet magne¬ 
tization and the electron injection parameters. 

There are two physical ingredients entering the mod¬ 
eling: 


• The first is the cooling of electrons and their syn¬ 


chrotron radiat ion, which is well understood (Uhm 


& Zhang 2014). This physical ingredient controls 


the shape of the spectrum. The general Band-like 
spectral-shape can be reproduced regardless of the 
details of the model parameters. 


• The second physical ingredient is related to the en¬ 
ergy dissipation physics and particle acceleration 
physics, which controls the model parameters such 
as Ho, 5, q , 7 i n j, and determine the Band parame¬ 
ters <a, E p , and /3. This likely invokes magnetic re¬ 
connection in a highly variable, highly magnetized 
environment, and the detailed physics is messy and 
poorly studiecQ 


As a result, one cannot predict from the first principles 
the values of those parameters, but can only infer them 
from the data. In any case, the inferred parameters as 
listed in Tabl e 1 fall i nto t he reasonable ranges from ex¬ 
pectation (Uhm & Zh ang| 2014), and the apparent corre- 


6 Recent numerical simulations (e.g.|Deng et al. 2015| are start¬ 
ing to tackle these difficult problems. 
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Figure 5. Comparison between the Band function fitting and synchrotron model fitting for Interval 7. The spectra have been re-binned 
into wider bins for clarity. Top two: observed count spectrum vs. model. Top Left shows the Band function fitting and Top Right shows 
synchrotron model fitting. Bottom two: the de-convolved photon spectrum plots. Bottom Left shows Band function fitting, red line is the 
theoretical photon spectrum of the Band function and data points are the “observed” photon flux which is obtained by de-convolving the 
observed count spectrum using instrument responses. Bottom Right: same as Bottom Left but for the synchrotron model fitting where the 
red line is the theoretical photon spectrum of the synchrotron model. 


lations among parameters are also generally understand¬ 
able. 


4.5. Non-existence of the Thermal Component 

For GRB 130606B, the Band function or the syn¬ 
chrotron model can adequately fit the data. There is 
no need to introduce a quasi-thermal component, which 
is presumably the emission from the photosphere of the 
outflow. The non-detection of the photosphere com¬ 


ponent poses a significant constraint on the GRB jet 
composition, which should ha ve a significant fraction of 
Poynting flux in the outflow (Zhang & Pe’er 2009 Gao 
fc Zhang||2015 ). 


5. SUMMARY AND DISCUSSION 


In this paper, we successfully fit the time-resolved spec¬ 
tra of GRB 130606B with a fast-cooling synchrotron radi¬ 
ation model that incorporates magnetic field decay. The 
same data are adequately fit by the empirical Band func- 





























E (keV) 


Figure 6. Comparison between the two theoretical photon models 
in Figure [5] 


tion. This is the first time that we give Band function a 
physical meaning through directly fitting the real GRB 
data with a physical model incorporating fast-cooling 
synchrotron radiation 


Our fast-cooling synchrotron model (Uhm & Zhang 


2014) involves a minimal number of physical assump- 


tions. It is impressive that the time-resolved spectra 
of GRB 130606B, which are characterized by the Band 
function with a wide range of spectral parameters, all 
can be equally well fit by the synchrotron model. An 
<a-index as hard as —0.52^0.24 can reproduced by the 
synchrotron model. 

Limited by the computational resources, it is beyond 
the scope of this paper to fit the spectra of a large sample 
of GRBs. With a 3.5 GHz 6-Core Intel Xeon E5 CPU 
on a Mac Pro workstation, it took about 160 CPU hours 
to finish synchrotron model fit for one time slice. Alto¬ 
gether it took about 4000 CPU hours to finish all the 
spectral fitting. The model fits to this particular GRB 
demonstrate that the synchrotron model can successfully 
explain typical Band function spectra and a systematic 
study of a much larger sample is highly encouraged for 
the future studies. 


We thank Xiaohong Zhao, Michael J. Burgess, Robert 
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helpful comments and suggestions. B.B.Z. acknowledges 
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acknowledges a NASA ADAP program (NNX14AF85G) 
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Figure 7. We compare the parameters of the physical model with those of the Band function and show the correlations between those 
pairs with strong and moderate correlations. The Pearson correlation coefficients of each pair are also marked in each plot. 














